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I Within the framework of three- flavor neutrino oscihations, we consider the 

\ physics potential of — > appearance and — > survival measurements 

at a neutrino factory for a leading oscillation scale ~ 3.5 x 10"'^ eV^. 
Event rates are evaluated versus baseline and stored muon energy, and opti- 
, mal values discussed. Over a sizeable region of oscillation parameter space, 

matter effects would enable the sign of (5m^ to be determined from a com- 
parison of i^e with event rates and energy distributions. It is 



\^ ' important, therefore, that both positive and negative muons can be stored in 

^ the ring. Measurements of the — > survival spectrum could determine 

the magnitude of dm^ and the leading oscillation amplitude with a precision 
of 0(l%-2%). 



I. INTRODUCTION 



The SuperKamiokande (SuperK) collaboration [Q] has published evidence that muon 
neutrinos, produced in the earth's atmosphere by cosmic rays, oscillate into other neutrino 
flavors. The z/^ survival probability in vacuum for neutrinos of energy (GeV) traversing 
a distance L (km) is given by: 

P{v^, Vf,) = sin^ 26'atm sm'^{l-2675mjrj^^L/ E^) . (1) 

The SuperK results for the oscillation amplitude sin^ 26'atm and the oscillation scale Sm\rj.-f^ 
are in accord with results obtained from other experiments 0. From the zenith angle 
distribution of the muon events, which is related to the distribution in L/E^, the mass- 
squared difference scale of the oscillations was inferred to be |jl]J^ 

(5miTM = (3.51^^) X 10-3 eV^ (2) 
and the amplitude of the oscillations was found to be maximal or nearly maximal, 

sin2 2^ATM = llo.2- (3) 

No zenith angle dependence was observed by SuperK for electron events, so it is con- 
cluded that electron-neutrinos do not undergo appreciable oscillations at the 5m\rj.y^ scale of 
Eq. (0) and that the muon-neutrinos oscillate dominantly to some neutrino flavor other than 
electron-neutrinos. This interpretation is consistent with the stringent lower limits on the 
electron-neutrino survival probability P{i'e '^e) from the CHOOZ reactor neutrino experi- 
ment . Thus it is inferred that atmospheric muon-neutrinos oscillate to tau-neutrinos or to 
a new sterile neutrino species that has no Standard Model interactions. In the latter case, 
matter effects would be expected to distort the zenith angle distributions at large zenith 
angles. The SuperK data presently disfavor muon-neutrino oscillations to sterile neutrinos 
by two standard deviations 

For Sm^ given by Eq. (|]), the flrst minimum in the survival probability occurs at 

L/E^ ~ 350tt?o km/GeV . (4) 

The detector averaged ^(t'^j v^) measurement by SuperK does not resolve this minimum 
because of the smearing over L and inferred Ey values. Accelerator-based experiments are 
thus essential to establish the existence of the minimum in the survival probability. The 
K2K 1^ experiment from KEK to SuperK, already underway, has L = 250 km and average 
neutrino energy (E^) = 1.4 GeV, with an L/E^ range of 125-250 km/GeV. The MINOS 
experiment from Fermilab to Soudan has a baseline L = 732 km and 3 beam options with 
{Ey) = 3, 6 and 12 GeV, giving L/E^ = 50-250 km/ GeV. Other experiments from CERN 
to detectors in Gran Sasso with a baseline L = 743 km are in the final design stages. The 
energy dependence of the charged-current rate, the neutral-current rate, and r-appearance 
with higher energy beams will test the neutrino oscillation hypothesis. The value of 
will be measured in these long-baseline experiments to higher precision than possible with 
atmospheric neutrinos. For example, MINOS may ultimately be able to make a 10% 
determination of 5m^. 
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Thus the next generation of accelerator long-basehne neutrino experiments are expected 
to firmly establish neutrino oscillations and improve our knowledge of Srn^. However, the 
dominant flux components in these accelerator experiments are and z/^. To make further 
progress in determining all of the parameters describing the oscillations it is desirable to 
also have z/g and Ve beams, in addition to beams of higher intensity. With this in mind, it 
has been pointed out that if an intense muon source of the type being developed for 
a possible future muon collider [|TU| is used together with a muon storage ring having long 



straight sections, the resulting "neutrino factory" would produce intense beams containing 
Vf. [Vf.) as well as {i-'fi)- In the original neutrino factory proposal [^] it was shown that if 
the storage ring was tilted down at a large angle, the neutrino beams would be sufficiently 
intense to produce thousands of interactions in a reasonable sized detector on the other side 
of the Earth. Neutrino factories have caught the attention of the community, and several 
groups are developing the concept with BNL, CERN, FNAL, and KEK all being considered 
as possible sites [|rT]-[T^. 



In Ref. [|14| we explored long-baseline neutrino physics at a neutrino factory within 
the framework of two-fiavor oscillations. In the present paper we extend our analysis of the 
physics potential at neutrino factories, presenting results within the framework of three-flavor 
oscillations. In particular we discuss: (i) the precise determination of ^Jtj-atm (~ 1% ~ 2%) 
from the measured i/^ — >■ survival probability; (ii) a proof that the sign of 5?Ti^rpj^ can in 
principle be extracted from v^. and z/^ v^. appearance measurements by exploiting 

matter effects that modify these oscillation probabilities; and (iii) the optimal LjE and 
methodology for the measurements. Table | summarizes the baselines, average electron 
density in the earth, and the dip angle for various possibilities. Our present analysis focuses 
on long-baseline experiments with L = 732 km (Fermilab — > Soudan), L = 2800 km and 
2900 km (Note: Fermilab ^ SLAG ~2900 km, Fermilab Seattle, Washington ~2700 km), 
and L = 7332 km (Fermilab Gran Sasso). 



II. LEADING NEUTRINO OSCILLATION 

There are experimental indications of neutrino oscillation effects from the LSND accel- 
erator experiment , from the atmospheric neutrino anomaly [l|J^,|l^ , and from the solar 
neutrino deficit |]T^|I8|. Three neutrino mass-squared differences are required to completely 



explain all these phenomena. However, three neutrinos provide only two distinct scales. 
Therefore, a sterile neutrino must be invoked if all the experimental indications are real. 
Since the significance of the LSND effect is not at the discovery level, a common approach is 
to set this anomaly aside until it is confirmed or rejected by the forthcoming Fermilab Mini- 
BooNE experiment and analyze the solar and atmospheric neutrino deficits within the 



framework of three-neutrino oscillations. This is the route followed in the present analysis. 
Thus we take the atmospheric sale to be the leading oscillation. 

With three neutrinos, the fiavor eigenstates z/q, (a = e, fi, r) are related to the mass 
eigenstates uj (j = 1, 2, 3) in vacuum by 

j 
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where U is the unitary Maki-Nakagawa-Sakata (MNS) mixing matrix [|20|. We parametrize 
Uhy 



U 



C13C12 

~C23Si2 — Si3S23Ci2e* 

1 iS 

\ ■523-512 ~ ■5l3C23Cl2e 



C13S12 



C23C12 ~ ■5l3'523'5l2e C13S23 
"■523012 — Si3C23Si2e C13C23 



(6) 



where cjk = cos 9 jk and Sjk = sin 6 jk- For Majorana neutrinos, U contains two further 
multiphcative phase factors that modify diagonal entries of U, but these do not enter in 
oscillation phenomena. 

The importance of forward scattering of neutrinos on electrons in the propagation of 
neutrinos through matter was first pointed out in Ref. [^. The existence of resonance 
effects in propagation of neutrinos through the earth in a constant density approximation 
was then discovered in Ref. p2]. Corrections to the magnitude and sign of the coherent 
amplitude given by Ref. were made in Refs. p2|,p3|. The effects of matter are prominent 



in neutrino oscillation solutions to the solar neutrino anomaly [T^jlMl ■ Early calculations 
for long-baseline experiments were presented in Ref. [^. Many long-baseline analytical and 
numerical studies have subsequently been made that include various refinements |26-p2|. 
The propagation of three neutrinos through matter is described by the evolution equation 



dx 



E 

13 



A ^ ^ 



2Ey j 2Ey 



1^(3 ■ 



(7) 



where x = ct and A/{2Ei,) is the amplitude for coherent forward charged- current scattering 
of z/p on electrons, 



A = 2V2GFYepE^ = 1.52 x IQ-^ eV^Y,p{g/cm^)E{GeV) 



(8) 



Here Ye{x) is the electron fraction and p{x) is the matter density. For neutrino trajectories 
that pass through the earth's crust, the average density is typically of order 3 gm/cm^ and 
Ye ~ 0.5. The propagation Eq. (0) can be re-expressed in terms of mass-squared differences 

as 



1 



1 = > 

dx y 2E^ 



Up. 



(9) 



This evolution equation can be solved numerically for given input values of the and 
mixing matrix elements. 

In our analysis we obtain numerical solutions of Eq. for three-neutrino oscillations 
taking into account the dependence of the density on depth using the Preliminary Reference 
Earth Model |Q . However, it is instructive to first consider an approximate analytic solution 
based on the (5m^ hierarchy 



|(5m32p ~ |5m3i|2 > |5m2ip. 



(10) 



implied by oscillation solutions to the atmospheric and solar neutrino anomalies. Then, for 
short enough distances, the Sm2i term in Eq. (|^) can be dropped and the other two (5m^ set 
equal. 
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(11) 



The evolution equations are then 



''dt 



( 



5rn? 
~2E 



V3 



a. 



t3\ 







1 









(12) 



The propagation matrix has eigenvalues 



Ai = 



6m'^ A 



6m? 



S 



5m? 

1e 



1 + 



A 
dm^ 



+ S 



(13) 



where 



S 



Y. 



A 



-cos 2^13 +sin2 2^ 



13 



(14) 



For propagation through matter of constant density, the flavor eigenstates are related to the 
mass eigenstates i/™ by 



where 



/ 



r" 



13 -^13 \ 

C23 ~'5i3'523 (^\-iS2Z 



(15) 



(16) 



■523 ~'5i3C23 C;^3C23 



and 6*^ is related to dxj, by 



tan 26* 



sin^ 2^13 
^'"cos2ei3-^ 



We note that has the form of the vacuum V with the substitutions 

^13^^r3, ^23^^23, ^12^7r/2, 5 = 0. 



(17) 



(18) 



There are only two mixing angles because the oscillations associated with the solar scale 
(5m2i) have not yet developed. The effective phase angle 5 vanishes, so CP violation is not 
possible in the leading oscillation, even after matter effects are included [^. Equation (|T^ 
implies that 



sin^ Id 



sin^ 2e 



13 



fA-cos2^^i3') +sin2 2ei3 



(^-cos26i3^ 
Thus there is a resonant enhancement for 

A = Sm^ cos 2^13 



(19) 



(20) 
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or equivalently 



E,,^15 GeV 



3.5 X 10-3 eV2 



1.5 g/cm^ 



COS 29 



13 



(21) 



The resonance occurs only for positive Sm"^. For negative the oscillation amplitude in 
(|T9|) is smaller than the vacuum oscillation amplitude. Thus the matter effects give us a way 
in principle to determine the sign of Sm^. 

For a constant Earth density profile, the transition probabilities for a given neutrino 
flavor are given by 



where 



Pi^a - = -4 E U^'.UZU^.UJl sin' A™. , 
j<k 



A™. = L{X, - A,)/2 



(22) 



(23) 



and the Xj are the eigenvalues of the neutrino matrix in Eq. ( |TB| ) . The transition probabilities 
in the leading oscillation approximation are 



sl. Sin' 2eZ sin' A™ , 
cl, sin' 2^r3 sin^ A™ , 



(24) 



sin^ 2^23 



sin'A^ + (sin C^)^ sin' A 



^21 



13; 



^31 



(sin^™ cos ^S)' sin' ^32 



Here the oscillation arguments are 

1 



Am 
^32 



Am 
^31 



Ao- 



1 + 



A 
5im? 



+ S 



21 



Ar 



1 + 



A 
5iv? 



-S 



where S is given by Eq. (|I1|) and 

Ao 



1.267 



5ml^ (eV2) L (km) 



(25) 



(26) 



AE " (GeV) 

The A^ term in P(z/^ — > Vr) must be retained here because it is not necessarily negligible 
compared to A^, due to matter effects. From Eqs. (0), ( pif ) and (pS]), we see in the limit 
^13 that ^'J^ — > and the Vf. transition probabilities vanish. We see from Eq. (|^) that 
matter effects are possible for z/^ — >■ Ur oscillations in a three neutrino model. However, 
these matter effects disappear in the ^13 = limit, where 



S 



A 



6m' 



ATM 



and 



P{^^ Vr) = sin' 2^23 sin' (Ao) 



(27) 



(28) 



To have substantial r event rates from z/^ i>r oscillations requires high energy neutrino 



beams because of the kinematic suppression of the r production cross section near threshold. 
The z/e — >• and z/^ oscillation probabilities can be studied with either low or high 

energy beams and we focus our attention on them in the rest of this paper. A non-zero value 
of is essential for the occurrence of these oscillations and for the observation of matter 
effects. 
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III. NEUTRINO FLUXES AND CROSS-SECTIONS 



The distribution n{x,Q) = (P'N/dxdVt of neutrinos from the decays of an ensemble of 
polarized negatively- charged muons in the muon rest frame is given by 

n,^ (x, ^) = ^ [(3 - 2x) + (1 - 2x)P^ cos 9] , (29) 

np,(x,fi) = ^[(l-x) + (l-x)P^cos^] , (30) 
47r 

where x = 2i?j,/m^, 9 is the angle between the neutrino momentum vector and the muon 
spin direction and is the average muon polarization along the beam direction. The 
corresponding distributions for and z/g from /z"*" — >■ e^v^j^Ve are obtained by the replacement 
P^- — >■ —Pfj,+ in Eqs. ( pQ]) and (|30D. Only neutrinos emitted in the extreme forward direction 
(cos^ ~ 1) are relevant to the neutrino flux for long-baseline experiments; in this limit 

= xE^ (31) 

in the lab frame. The flux at a distance L from the storage ring can be approximated by 

(32) 

where 7 = E^/m^ and no is the number of neutrinos (or antineutrinos) in a given beam. The 
charged- current (CC) interaction cross-sections in the detector grow approximately linearly 
with the neutrino energy 

a^M ^ 0.67 X 10-3^ cm^ x E^, (GeV) . (33) 
apTv ^ 0.34 X 10"^^ cm^ x (GeV) . (34) 

Thus, the event rates in the absence of oscillations are 

N^{E^f/L\ (35) 

The z/^ ip^) event energies peak at the stored /i~ (/i"*") energy while the Pe (t'e) energies peak 
at about 3/4 of the stored /i~ (/i"*") energy. The lepton spectra from the charged-current 
interactions can be obtained by folding the event rates with the da/dy distributions, where 
y = l- Ee/E,. 

For no oscillations, the average observed neutrino energies would be 



=0.7E,-, (36) 
{E,J = O.QE^- . (37) 

Thus, from Eq. (|^), the first maximum in the leading oscillation occurs for 

L/E^ ~ 230l?J|] km/GeV . (38) 

For maximal oscillation effects the baseline and stored energy should be chosen to satisfy 
Eq. (psp. In order to have sufficient event rates at > 1000 km distances, the lowest practical 
muon energy is about 10 GeV. The corresponding L from Eq. (|3^) is close to the 2900 km 
baseline from Fermilab to SLAG. 
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IV. OSCILLATION CHANNELS 



The decays of stored muons along the straight sections of the storage ring give well- 
coUimated neutrino beams of specified flavor composition. For stored the fi~ —>■ t'^z/eC" 
decays give and z/g beams with the known energy distributions given by Eq. (^) (see 
Fig. 4 of [Q)- In this paper we restrict our attention to detection of muons. The relevant 
oscillation signatures for stored /x~ are given in Table ||. Corresponding signatures for stored 
/i"*" are obtained by the interchange u ^ v and changing the signs of the electric charges of 
the leptons. 

The oscillations to Vj. feed the same final state as oscillations to z/^ through the decays 
r ^ /i of the produced r-leptons. Thus detection of muons determines linear combinations 
of oscillation probabilities 



N{^x-) = (<l>P(z/^ ^ z/^)a(z/^ ^^^')) + {^P{i^^. yr)o{vr ^ t')) BF(r- ^ /x") , (39) 
iV(/x+) = (<l>P(z7e ^ z/,Ja(z7^ ^ /"+)) + L'^Pipe ^>r)a(z/, ^ r+)) BF(r+ ^ /i+) , (40) 



where $ is the flux and ( ) denotes a spectrum average and integration over final state 
energies. Neglecting the small matter effects in z/^ — z/^, the leading oscillation rates from 
Eqs. dl) and (|8D would be 

Ar(/i-) = ($ [l - sin^ 2^™ sin^ A!^2 - sin^ 2^23 sin^ Aq 



+ sin' 2^23 sin' A32 o{y, ^ r") ) BF(r- ^ /x') , (41) 
iV(/i+) = sin' 2^r3 \^ sin' A^2 ^(^,. ^ Z^' 

,2 o/am / /IS „;„2 



+ sin' 2^^ ($ sin' A'3^2 ^('^r ^ r+)) BF(r+ ^ ^+) . (42) 

The r-contributions will be suppressed at low stored muon energies (~10 GeV) by the kine- 
matic suppression of the tau cross section near the threshold for tau production. The tau 
neutrino interactions can be identified by (i) kinematics (p^ relative to the neutrino beam 
direction and the ^/-distribution of the decay muon), or (ii) direct evidence for r 1-prong 
or r — s> 3-prong decays by detecting kinks or displaced vertices in, for example, emulsion 



detectors, or by imaging Cherenkov rings that the tau generates in Cg F14 liquid |35]. Hence 



forth we concentrate our analysis strictly on primary muons from z/^ under the assumption 
that the contributions of muons from tau decay can be resolved. 



V. PREDICTED NEUTRINO OSCILLATION EVENT RATES 

For a quantitative analysis of neutrino oscillation event rates in long-baseline experi- 
ments, we begin with the following set of oscillation parameters: 





= 3.5 X 10" 


-3eV' 




|5m'i| 


= 5 X 10"^ 


eV' 




Sl3 


= 0.10 


(sin' 2^13 


= 0.04) 


■S23 


= 0.71 


(sin' 2^23 


= 1.00) 


S12 


= 0.53 


(sin' 2^12 


= 0.80) 


5 


= 
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The values of ^23 govern the atmospheric neutrino oscillations. The values of Sm^i 



and S12 control the solar neutrino oscillations. We select the large-angle matter solution ||36 
because it is the most optimistic choice for obtaining any effects in long-baseline experiments 
from the subleading Sm^i scale. The value of S13 determines the size of the matter effects 
for the leading oscillations. The value of Sm^i is determined from 6771^2 and 6mli by the 
sum rule 

Smli = 5m\2 + 5m\^ . (44) 



The subleading oscillation effects associated with the 5m2i in Eq. (|4^) are usually small in 
comparison with the leading oscillation, so either of the other two solar solution possibilities 
(small angle matter or large-angle vacuum oscillations) would lead to results similar to the 
representative case above; we will remark below when the subleading oscillation effects are 
appreciable. The long-baseline oscillations are controlled by the three parameters 5m'^2-i "523 
and S13. The matter effects in the leading oscillation are critically dependent on S13. 

We calculate neutrino event rates, with and without oscillations, for a neutrino factory 
with 2 X 10^° stored muons/year and a 10 kiloton detector. For the neutrino flux, energy 
distributions and cross sections we use Eqs. (p9D-(p4D. The oscillation probabilities are 
calculated by integrating Eq. (^ numerically from the source to the detector using a Runge- 
Kutta method, and averaging over the neutrino energy distribution. For the Earth's density 
we use the results of the Preliminary Reference Earth Model |^3[. In these initial exploratory 



results we neglect the stored muon beam momentum-spread and angular-divergence, which 
we later incorporate in Sec. Fig. |l|a shows the event rates at L = 732 km versus stored 



muon (or antimuon) energy for the oscillation parameters in Eq. (|43|) . We show in Fig. |T]b 
the corresponding event rates with dm^-^ ~ 10~^° eV^ for the subleading oscillation scale, 
i.e., appropriate for vacuum long- wavelength solar oscillations. We note the following: 

(i) For the appearance channels u^, — > z/^ and z/g — >■ z/^, the main difference in rates 
comes from the difference in the neutrino and antineutrino cross sections. Matter 
effects, which arise when the sign of 5m\2 is changed, are seen to be relatively small in 
Fig. |l|b. This is to be expected since the characteristic wavelength for matter effects 
at these energies is of order the Earth's diameter, which is much larger than 732 km. 

(ii) The effects of the subleading mass scale, 6m2i, are evident by comparing the ap- 
pearance channels in Fig. |l]a to those in Fig. |I|b. For the large-angle MSW solution, 
Sm2i = 5 X 10~^ eV^, the interplay of the subleading scale with matter and the sign 
of 5^32 affects the rates at the 20% level. We have found that this effect increases as 
sin^ 2^13 0, when the subleading leading oscillation, with its larger oscillation ampli- 
tude, can begin to compete with the leading oscillation, which has a small amplitude 



in this limit p7 |. 



(iii) At lower muon energies (10 GeV or less) the survival probabilities for the z/^ —>■ and 
i^fj. i^fj. channels decrease with decreasing energy, but do not reach their minimum 
values, on average, for > 5 GeV. 

(iv) The survival rates do not depend appreciably on the sign of ^mgg, which means that 
matter effects are small in these channels at this distance, as is the case for the ap- 
pearance channels. 
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In Fig. ^ we show similar results for L = 2900 km, from which we conclude: 

(i) The largest event rate suppression in the survival channels occurs at energies E^^ ~ 
10 GeV. 

(ii) Changing the sign of causes noticeable changes in the event rates in the appear- 
ance channels and z/g Ufj,. By comparing Fig. ^ with Fig. |^ we see that 
the subleading oscillation scale has a smaller relative effect for L = 2900 km than for 
L = 732 km (Fig. |T]), so these changes must largely arise from matter effects. 

In Fig. ^ we show similar results for L = 7332 km, from which we conclude: 

(i) Matter effects have now grown quite large, as evidenced by the large variation in the 
appearance event rates in Fig. |^ when the sign of (5m|2 is changed. We have found 
that the results using the large-angle MSW and vacuum solar scales for the subleading 
oscillations are nearly identical at this distance. 

(ii) A minimum in the survival probability is clearly visible near = 25 GeV. Also, the 
survival rates are now sensitive to the sign of -^'^i ~ 15 — 25 GeV, indicating 
matter effects. Unfortunately, the number of events in one of the appearance channels 
(either z/g Ufj, or z/g — > z/^, depending on the sign of f^^ls below 10 for these 
energies, even for sin^ 2^13 = 0.04. In general, the event rates at 7332 km are below 
those at 2900 km due to a lower flux. 

(iii) By comparing Fig. ^a with Fig. ^b we see that matter effects in the survival channels 
are drastically reduced as sin^ 26^13 is decreased, as predicted in Eq. (^81). 

(iv) The number of events in the appearance channels is sufficiently large to allow a deter- 
mination of the sign of 5m|2, at least for the L = 2900 km and 7332 km distances. 

Next we show in Fig. ^ the L dependence of the oscillations at a fixed muon energy of 
10 GeV. Here we see that L ~ 2900 km is the first minimum in the survival probability and 
that L ~ 5500 km corresponds to a maximum in the survival probability; matter effects are 
sizable for L ~ 5000 — 8000 km. In the appearance channels, matter effects turn on as L 
increases, causing one of the appearance channels to be highly suppressed. 

Similar results for E^^ = 50 GeV are shown in Fig. We see that the appearance 
event rates scale roughly with E^, as expected from the product of the E^ behavior of the 
unoscillated rates (Eq. (|55|)) and the E~'^ dependence of the oscillatory factor for small 
oscillation arguments (Eqs. p^-p^). The L dependence is relatively flat at low L for the 
channel that is not suppressed by matter, as expected from the product of the 1/L^ behavior 
of the flux and the dependence of the oscillatory factor. The survival rates are much 
higher for 50 GeV than for 10 GeV, but they do not reach a minimum for L < 10000 km. 

Figures ^|-^ assume I^Ti^lgl = 3.5 x 10~^ eV^, the favored value from atmospheric mea- 
surements from SuperK. Long-baseline experiments should be able to determine l^rrigjl with 
high precision if the event rates associated with oscillations are appreciable. Figure ^ shows 
the sensitivity of the event rates to changes in 5m|2 for L = 2900 km and = 10 or 50 GeV. 
At En = 10 GeV the survival probability has a minimum near l^mlgl = 3.5 x 10~^ eV^, but 
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the appearance rates are modest. At i?^ = 50 GeV the appearance rates are much higher, 
but the survival probabihty has a minimum only for higher | ^mgg | . Therefore we see that at 
L = 2900 km, if 5m|2 is close to the value favored by SuperK, there is a trade-off between 
choosing a relatively low to optimize the survival probability measurements, or choosing 
a higher energy to optimize the appearance signal. 

Figure ^ shows similar curves for L = 7332 km with i?^ = 25 and 50 GeV. At 
Efj, = 25 GeV the survival probability has a minimum near l^m^jl = 3.5 x 10"'^ eV^, 
with appreciable appearance rates for nearby values of 5m|2, and matter effects are evident 
for |(5ni32| > 3 X 10~^ eV^. At = 50 GeV, the appearance rates are higher than for 
25 GeV (especially in the matter-suppressed channel) and the disappearance probability is 
appreciable over a range of l^n^lgl- However, matter effects in the survival probability are 
noticeable only for I^Ti^lgl > 7 x 10"'^ eV^. 

Figure |^ shows the minimum value of sin^ 26'i3 that gives 10 z/g — z/^ appearance events 
versus E^ for L = 732, 2900, and 7332 km, when Sm^i ~ 10~^° eV^ (the solar vacuum 
oscillation value). The lower (upper) curves are for matter-enhanced (matter-suppressed) 
oscillations with 6ml2 > (5m|2 < 0); the corresponding curves for z/g —>■ appearance 
can be empirically obtained by changing the sign of 5m|2 and multiplying by 2. We see 
that increased E^ improves the ability to discover appearance channels, although there are 
limiting returns for going to very high E^. The sensitivity to the matter-suppressed channel 
is especially weak at longer distances (e.g. L = 7332 km in this case). The corresponding 
curves for the large-angle MSW solar solution (e.g. Sm^i = 5 x 10~^ eV^) are much different 
for L = 732 km, and somewhat different at L = 2900 km, due to contamination by the 
subleading oscillation, as expected from the discussion of Figs. |1] and |^. If large-angle solar 
oscillations are the correct solution to the solar anomaly, then L = 7332 km has the merit 
that the subleading oscillation does not affect the interpretation of z/g z/^ appearance, as 
it does at L = 732 and 2900 km. 

As discussed in Sec. 0, the leading oscillation approximation leads to a simple expressions 
for the oscillation argument of the z/g — >■ z/^ probability (see Eqs. ( pSj) and (pT])). Then the 
^^J. appearance event rate can be approximated by 



N ~ (<l>P(z/e z/^)a(z/^ ^ fi-)) (45) 




where the angle brackets denote averages over neutrino energy, and (E^) = 0.6i?^. The 
average probability can be estimated by evaluating the probability at (E,^) in the small 
sin^ 29i3 limit. 



sin^ 26*13 . 2 



1 _ (^) ^ 



sm 



[l_27^^32^ 


1 - 


{A) 1 









(46) 



where (A) is given by Eq. (H) evaluated at {Eu). Equation (^) gives reasonably accurate 
results whenever sin^ 2^13 is below 0.01 provided that contributions to the subleading oscilla- 
tion are not important. Approximate antineutrino rates can be obtained by the substitution 
A —A in Eq. ( ^61) and by dividing the results of Eq. ( ^5)) by 2. 
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VI. MORE DETAILED SIMULATIONS 



To obtain a more realistic calculation of the event rates at a neutrino factory we have 
simulated a neutrino beam formed by muons decaying along a 1 km long straight-section. 
In our simulations the muon beam has a momentum spread given by (Jp/p = 0.02, and 
horizontal and vertical beam divergences given by = O.I/7 where 7 = E^/m^. The finite 
beam momentum spread and angular divergence modify the event rates at a distant site by 
a few percent. We have also included in our simulations a parametrization of the detector 
resolution functions for an iron-scintillator neutrino detector @J^. We assume that muon 
energies can be measured by range with a precision given hy (Te/ E = 0.05. Shower energies 
are assumed to be measured with precisions given by cte/E = 0.53/-\/E if > 3 GeV and 
aE/E = 0.80//E if E < 3 GeV with, in both constant term of 0.07 added in 

quadrature. 

For each of the relevant oscillation channels, the predicted annual CC event samples (from 
2 X 10^° muon decays) corresponding to the oscillation parameters of Eq. (^) are shown in 
Table |Tl| for both signs of Sm'^^- The tabulated rates are for a 10 kt detector at baselines 
of L = 732 km, 2800 km, and 7332 km, with = 10, 30 and 50 GeV. The statistical 
uncertainties on the calculated event rates from the Monte Carlo simulations are about 2.5%. 
Consider first the sensitivity of the oscillation signals to matter effects, and hence to the 
sign of 5ml2- Note that most of the event rates are insensitive to matter effects. However, 
the z/g — > CC rates are significantly modified by matter, and hence a measurement of 
events with wrong-sign muons can in principle determine the sign of 6m'^2- To illustrate 
this, consider the Ue — > z/^ CC rates when L = 2800 km and E^ = 50 GeV (30 GeV). If 
6m'^2 is positive, 894 (486) wrong-sign muon events are expected in a 10 kt-yr data sample. 
However, if 6m'^2 is negative, the corresponding rates are 338 (130) events. Hence, if all the 
oscillation parameters are known except the sign of 6771^2, the event rates can be used to 
determine the sign. The statistical significance of this sign determination improves slowly 
with increasing Efj^. In practice, the oscillation parameters will not be precisely known. 
However, the sign of 5m|2 can still be uniquely determined from the wrong-sign muon 
measurements if both z/g z/^ rates (positive stored muons) and Ue i^ij. rates (negative 
stored muons) are determined. For the examples we are considering, if ^rrigg is positive the 
i^e —>■ i^fi rates will be larger than the z/g — > z/^ rates by a factor of a few. In contrast, if Sml2 
is negative the z/g —>■ rates will be comparable to, or smaller than, the z/g — > z/^ rates. 
There is additional information in the spectrum of CC events producing wrong-sign muons. 
The predicted measured energy distributions (including detector resolutions), for CC events 
containing wrong-sign muons are shown as a function of both the magnitude and sign of 
5m\2 in Figs. ^ and ^ when, respectively, 30 GeV positive and negative muons are stored 
in the ring, and L = 2800 km. Both the shapes and normalizations of the distributions are 
sensitive to the sign and magnitude of 5m|2- In particular if 5m|2 is negative (positive) the 
fi~ appearance events from z/g z/^ oscillations will peak at higher (lower) energies than 
the /i^ appearance events from z/g z/^ oscillations. We conclude that if both positive 
and negative muons can be stored at different times in the neutrino factory, then wrong- 
sign muon appearance measurements can distinguish the sign of (Jm^g and determine its 
magnitude. 

We next consider a less favorable region of parameter space in which the z/g 
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oscillation amplitude is reduced by a large factor. Table |^ presents the event rates when 
sin^ 26^13 is reduced by an order of magnitude (sis = 0.032 or sin^ 26*13 = 0.004) with the 
other oscillation parameters as in Eq. (^). Returning to our example (L = 2800 km, 
= 50 GeV) we note that if ^mgg is positive (negative) we expect 108 (25) wrong-sign CC 
events per 10 kt-yr when positive muons are stored, and 29 (30) events when negative muons 
are stored. Hence, in this case we are still above, but close to, the threshold below which 
a statistically significant determination of the sign of 61^1^2 ^il^ only be possible with larger 
data samples. Based on Tables |T| and |^ we can estimate the minimum values of sin^ 26*13 
(when all of the other parameters are as specified in Eq. (^)) for which (i) we expect a 
measurable wrong-sign muon signal (10 event sensitivity), and (ii) from the measured ratio 
of wrong-sign muon rates obtained when positive and negative muons are alternately stored 
we can determine the sign of {^^ 3 standard deviations). Figure |10| summarizes our 

estimates for these minimum sin^ 2^13 versus L and E^. Note that to measure a wrong-sign 
muon signal over the largest sin^ 2^13 range, high E^ and "short" L are preferred. However, 
if L is too short the unoscillated CC event rates get very large, and hence backgrounds 
become significant. If the background rates are at the level of 10~^-10~^ of the total CC 



rates ||3^ , then L = 2800 km is probably preferred over 732 km. To be able to determine the 
sign of o'^^r the largest sin^ 2^13 range, high E^ is once again preferred. Since matter 
effects become small for short L and event statistics become small for very large L, there 
appears to be an optimal L for determining the sign of Sm'^2- Of the 3 baselines we have 
considered, L = 2800 km is preferred. From these considerations we would conclude that 
for wrong-sign muon measurements L = 2800 km is a good choice for the baseline, and high 
energy stored muons {E^ = 50 GeV) are preferred, although we note that decreasing i?^ 
from 50 GeV to 30 GeV decreases the sensitive sin^ 2^13 range by less than a factor of two. 

We next turn our attention to muon survival measurements. The z/^ rates shown 

in Table |ITI| are significantly less than the corresponding rates in the absence of oscillations. 
As an example, consider L = 2800 km with E^ = 10 GeV and 30 GeV. The predicted mea- 
sured energy distributions (including detector resolutions) for interacting muon neutrinos 
are shown for the two energies in Figs. 11 and 12 respectively versus (5m|2 with the other 
oscillation parameters given by Eq. (^3]). The shapes and normalizations for the predicted 



distributions are very different from the expected distribution that would be observed in 
the absence of oscillations. The shapes of the predicted z/^ CC interaction distributions are 
sensitive to the magnitude of Sm'^2- The dip in the predicted rate that corresponds to the 
first oscillation maximum (when 1.267 5ml2L/E^ = tt) can be clearly seen. At the higher E"^ 
the contribution from poorly measured higher energy z/^ CC interactions reduces the signif- 
icance of the dip. To understand the statistical precision with which 6171^2 sin^ 2^13 can 
be extracted from a fit to the measured z/^ CC interaction distribution, we have generated 
Monte Carlo data sets with the appropriate statistics for various values of the oscillation 
parameters, and fit the resulting simulated distributions. As an example. Figs. |T3| and ^ 
show fit results for E^ = 30 GeV with L = 2800 km and E^ = 50 GeV with L = 7332 km. 
The precision of the fits for various L and i?^ are summarized in Table |V] for sin^ 26*23 = 1-0 
and Sm^2 = 0-35 x 10^ eV^. In order to extract the fitted information, we have had to limit 
the fits to energy ranges where the information is maximally available. It can be seen from 
Table that the parameters are ill constrained for the shorter baseline of 732 km. For the 
10 GeV ring, L = 2800 km and for the 30 or 50 GeV ring L = 7332 km seem to be the 
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optimal baseline lengths for maximum precision in sin^2^23 and 5m 



32- 



VII. CONCLUSIONS 

Within the framework of three-flavor oscillations, we have explored the sensitivity of 
muon appearance and muon disappearance measurements at a neutrino factory in which 
there are 2 x 10^° muon decays per year in a beam-forming straight section that points at 
a 10 kt detector. Three stored muon energies (10, 30, and 50 GeV) and three baselines 
(732, 2800, and 7332 km) have been explicitly considered. Our results are summarized in 
Table VI. 

If data can be taken alternately with positive and negative muons stored in the ring, a 
measurement of the wrong-sign muon appearance rates and spectra can uniquely determine 
the sign of 6m\2 provided the oscillation amplitude is sufficiently large (sin^20i3 > 0.005 for 
L = 2800 km and = 30 GeV). To demonstrate this, we have considered determining 
the sign of 5m|2 by measuring the ratio of wrong-sign muon rates when alternately positive 
and negative muons are stored. Of the three baselines we have considered, L = 2800 km 
is preferred for this measurement since L = 732 km is too short (matter effects too small) 
and L = 7332 km is too long (statistics too limited) to obtain good sensitivity to the sign of 
6m'^2 ■ We note that a global fit to both wrong-sign muon event rates and energy distributions 
with 6m'^2 (including its sign) left as a free parameter would improve the sensitivity to the 
sign of 5m|2, and may result in a different preferred baseline choice ||3^. The sensitivity to 
wrong-sign muon appearance, characterized by the minimum sin^2^i3 for which a signal can 
be detected at the 10 events per year level, improves linearly with E^^. For L = 2800 km and 
Efj^ = 30 GeV the minimum sin^26'i3 = 0.0007. This can be improved to sin^2^i3 = 0.0004 
by either increasing E^ to 50 GeV, or by increasing the detector mass by a factor of 1.8. 

The sensitivity for the survival measurements can be optimized by choosing L and -E^ 
so that L/ Ey is close to the first minimum in the survival probability. For L = 2800 km this 
corresponds to i?^ ~ 10 GeV. For 6771^2 = 0.0035 eV^/c^ and sin^2^23 = 1 the statistical pre- 
cision of the resulting 5m|2 and sin^26'23 measurements based on fitting the observed z/^ z/^ 
spectra are respectively 2.4% and 1.1%. However, optimization for the appearance channels 
suggests choosing a higher E^j,. For example, choosing E^ = 30 GeV the precisions for the 
6m'^2 and sin^26'23 measurements become 3.2% and 2.0%. With these levels of statistical pre- 
cision, systematic effects (for example, the uncertainty on the neutrino flux) may dominate. 
For the region of three-flavor-mixing parameter space that we have explored, we conclude 
that L = 2800 km with E^j^ = 30 GeV would enable the very precise determination of 1 5^^321 
and sin^ 2^23 from muon survival, the determination of sin^ 2^13 from muon appearance, and 



My' 



the sign of Sm'^2 from matter effects (e.g., by comparing z/g with z/g — >■ z>^ 

The above analysis assumes that CP is conserved, i.e., 5 = 0. CP violation may be 
important whenever the effects of the subleading mass scale 5m2i are appreciable, e.g., at 
short and intermediate distances A detailed analysis of the CP- violating case will be 



given elsewhere [39 



Finally, we note that if no appearance signal were observed, the implied very low value 
of sin^26'i3 might provide a window of opportunity for detecting oscillations driven by the 
smaller (solar neutrino deflcit) scale Sm^i, should it be the large-angle MSW solution. This 
requires further study. 
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TABLES 



TABLE L Some candidate neutrino factory accelerator sites and long-baseline detector sites. 



Sources 



FNAL 




BNL 






CERN 




KEK 


Targets L (km) (Y^p) .^^^^H 


L(km) 




angle ^ ^ 


L (km) 


^ep) ariglc(°) 


L (km) 


a''n|ie(°) 


MINOS 732 1.67 3.3 


1715 


1.67 


7.7 


6596 


2.01 31 


8559 


2.20 42 


BNL 1289 1.67 5.8 








5921 


1.91 28 


9607 


2.31 49 


SLAC 2913 1.67 13 


4513 


1.69 


19 


8590 


2.20 42 


7720 


2.12 37 


Gran Sasso 7332 2.09 35 


6565 


2.00 


31 


743 


1.67 3.3 


8844 


2.23 44 


rvamioKa yioo z.zo 4o 


9667 


2.31 


49 


8774 


2.22 44 


252 


1.67 1.1 


TABLE IL Signatures for 


various 


oscillation channels when negative muons are stored. 












produced 




detected 


signature 


oscillation 






lepton 




muon 


z^jit-survival 
















z^r-appearance 
















i/^-appearance 
















^T-appearance 
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TABLE III. Predicted event rates in a 10 kt detector per 2 x lO^'^ muon decays, with the 
oscillation parameters specified by Eq. (|43|). 





sign of 




Ep+ (GeV) 




E^,- (GeV) 
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^n 


10 


30 


50 








r 

/ — 


■ 1 oz Km 
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qon nnn 
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i^fi t^fi 


+ 


11,200 
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1 z,uuu 
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22,800 


750,000 


3,520,000 




— 


11,200 


'i7'y nnn 


1 7Rn nnn 


22,900 


751,000 


3,530,000 


Ve (no osc) 




24,200 


Ri^fi nnn 

000, uuu 


'i n^n nnn 


12,400 


329,000 


1,520,000 


l^e J^fM 


+ 


158 


404 


1 000 


69.4 


276 


486 




— 


106 


41 n 


71 S 


65.2 


223 


380 








r 

— 


9snn Vm 
zouu Km 








(no osc) 




1140 


97 f^nn 
Z ( ,ouu 


1 9^ nnn 
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2200 


52,500 


244,000 


J^fM v^l 


+ 


168 


1 s snn 

io,oUU 


1 n7 nnn 
lU / ,uuu 


274 


36,000 


212,000 




— 


164 


1 s onn 


1 n7 nnn 
lU / ,uuu 
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36,200 


212,000 


fe (no osc) 




1400 


A'i '^nn 


91 1 nnn 

Zil,UUU 


900 


22,600 


102,000 


l^e J^fi 


+ 


54 


486 


894 


4 


88 


224 




- 


3.4 


130 


338 


36 


214 


356 








L = 


7332 km 








(no osc) 




152 


3880 


17,500 


294 


7660 


35,600 




+ 


84 


592 


7220 


130 


1040 


14,900 






72 


558 


7340 
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1130 


15,000 


Ue (no osc) 




21.8 


6680 


30,500 


120 


3300 


15,200 


l^e J^fi 


+ 


21.6 


274 


300 


0.66 


9 


19.8 






1.14 


18 


38 


12.8 


132 


145 
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TABLE IV. Predicted event rates in a 10 kt detector per 2 x 10^° muon decays, with 
sin^ 2^13 = 0.004 and the other oscihation parameters specified by Eq. (]43|). 
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TABLE V. 


Statistical 


precision 


in determining sin^ 2^23 and 


if sin^ 26*23 


= 1.0 and 


(5m32 = 3.5 X 


10-3 eV^. 


Fit result 


s are tabulated for various 


muon storage ring energies and 


basehnes. 
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TABLE VI. Summary of sensitivity versus baseline and stored muon energy. 
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FIGURES 
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FIG. 1. Muon survival and appearance events per 10 kt-yr versus muon energy at L = 732 km, 
assuming 2 x 1020 muon decays per year. The oscillation parameters are those given in Eq. (^) and 
(a) dm2i = 5 X 10~^ eV^ and (b) 5m2i 1 x 10~^^ eV^. The upper (lower) solid curve shows the 
rates without oscillations for neutrinos (antineutrinos). For the rates with oscillations, the results 
are shown for both signs of 6m'^2- 
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FIG. 2. Same as Fig. || except L = 2900 km. 
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FIG. 3. Same as Fig. except L = 7332 km, bm\ = 5 x lO'^ eV^, and (a) sin^ 2013 = 0.04 
and (b) sin^ 26*13 = 0.004. 
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FIG. 4. Muon survival and appearance rates per 10 kt-yr versus detector distance for (a) 
= 10 GeV and (b) = 50 GeV. The oscillation parameters used are given in Eq. 
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FIG. 5. Muon survival and appearance rates per 10 kt-yr versus \dml2\ at L = 2900 km for 
Ef^ = 10 GeV and (b) = 50 GeV. The other oscihation parameters are given in Eq. (|43|). 
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FIG. 6. Same as Fig. || except L = 7332 km. 
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FIG. 7. Minimum value of sin^ 2^13 that gives 10 Ve i^^i appearance events in a 10 kt 
detector with a source of 2 x 10^'^ muon decays/yr, for L = 732 km (sohd curves), 2900 km (dashed 
curves), and 7332 km (dotted curves). The oscihation parameters are given in Eq. ([43|) and 
1 X 10"^" eV^. The lower (upper) curve in each case is the result for 6171^2 > {^^32 < 0)- 



The corresponding curves for 



appearance are given approximately by interchanging the 



upper and lower curves and dividing by a factor of 2. 
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FIG. 8. Predicted measured energy distributions for CC events tagged by a wrong-sign (neg- 



ative) muon from Vf. 



oscillations, shown for various bm%^2i ^ labelled. The predictions 



correspond to 2 x 10 decays, 
and 8 given in Eq. (|4F 
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FIG. 9. Predicted measured energy distributions for CC events tagged by a wrong-sign (pos- 
itive) muon from — > i)^ oscillations, shown for various (^rrigj, as labelled. The predictions 
correspond to 2 x lO^'^ decays, = 30 GeV, L = 2800 km, with the values for 5mf2, S13, 523, •Si2) 
and 6 given in Eq. (|4^). 
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FIG. 10. The value of sin^ 26*13 that yields, in a 10 kt detector at L = 2800 km, (a) 10 events 
per 2 X 10^"//"*" decays (boxes), and (b) a three standard deviation determination of the sign of 5m'^2 
(circles) based on the ratio of wrong-sign muon rates when alternately 2 x 10^'^ /x^ and 2 x 10^*^ 
/i~ decay in the neutrino beam-forming straight section. The sin^ 2^13 sensitivity is shown versus 

and L (as labelled). The calculations assume the values for |(^?Ti32|, (^?Tif2, S23, S12, ^ given in 
Eq. (^3|). The shaded region is excluded by the existing data. 
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FIG. 11. Predicted measured energy distributions for CC events shown for four 

different 5m'^2 (darkly shaded distributions) as labelled. The predictions correspond to 2 x 10^*^ 
decays, E^j^ = 10 GeV, L = 2800 km, with the values for 5ml2, sn, S23, si2j and 5 given in Eq. (^3|). 
The predicted distribution has been used to generate a Monte Carlo dataset with the statistics 
corresponding to a 10 kt-yr dataset (points with error bars). The lightly shaded histograms show 
the predicted distributions in the absence of oscillations. 




FIG. 12. Predicted measured energy distributions for CC events shown for four 

different (5m|2 (darkly shaded distributions) as labelled. The predictions correspond to 2 x 10^*^ 
decays, = 30 GeV, L = 2800 km, with the values for (5m^2) ■^iS; ^23, 512, and 6 given in Eq. (^3|). 
The predicted distribution has been used to generate a Monte Carlo dataset with the statistics 
corresponding to a 10 kt-yr dataset (points with error bars). The lightly shaded histograms show 
the predicted distributions in the absence of oscillations. 
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FIG. 13. Fit to muon neutrino survival distribution for = 30 GeV and L = 2800 km for 
10 pairs of sin^20, dm? values. For each fit, the Icr, 2cr and Scr contours are shown. The generated 
points are indicated by the dark rectangles and the fitted values by stars. The SuperK 68%, 90%, 
and 95% confidence levels are superimposed. Each point is labelled by the predicted number of 
signal events for that point. 
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FIG. 14. Fit to muon neutrino survival distribution for = 50 GeV and L = 7332 km for 
10 pairs of sin^ 26, 5m^ values. For each fit, the la, 2a and 3a contours are shown. The generated 
points are indicated by the dark rectangles and the fitted values by stars. The SuperK 68%, 90%, 
and 95% confidence levels are superimposed. Each point is labelled by the predicted number of 
signal events for that point. 
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